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Abstract. We consider the thermodynamic behavior of a disordered interacting electron system in two 
dimensions. We show that the corrections to the thermodynamic potential in the weakly localized regime 
give rise to a non monotonic behavior of the specific heat both in temperature and magnetic held. From 
this effect we predict the appearance of adiabatic hysteresis in the magnetoconductance. Our results can 
be interpreted as precursor effect of formation of local moments in disordered electron systems. We also 
comment on the relevance of our analysis in three dimensional systems. 
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The recent discovery of a B = metal-insulator tran- 
sition (MIT) in high mobility silicon MOSFETs has re- 
vived the interest in disordered electron systems. Whereas 
conventional scaling theory predicts no metallic state in 
two dimensions Q the experiments show a metallic tem- 
perature dependent resistivity with a characteristic expo- 
nential behavior. 

Although there have been attempts to interpret these 
experimental findings within a phenomenological single 
parameter scaling theoryj^j, there is growing experimen- 
tal evidence as witnessed by magnetic field measurements 
jl| that electron-electron interaction in the spin channel 
is relevant and needs explicit consideration. Based on the 
existing theory of disorder and interaction^), it has been 
pointed out recently j(| that the very existence of a metal- 
lic state at low temperatures is indeed possible because 
of enhanced spin fluctuations, though a more thorough 
understanding of the behavior of various physical quanti- 
ties is doubtless required. In particular it has been sug- 
gested that magnetoconductance in parallel field and tun- 
neling measurements should provide good diagnostic tests 
for this theory. 

A further distinctive prediction of the theory of dis- 
order and interaction concerns the anomalous low tem- 
perature behavior of the thermodynamical quantities @,||. 
This appears to be consistent, at least qualitatively, with a 
number of experimental results, among which NMR, spin 
susceptibility, and specific heat measurements in 3d Si:P 
metallic samples^. At phenomenological level, a two fluid 
picture consisting of a free electron gas coupled to a set 
of localized magnetic moments seems to capture the main 
features of the experimental results. 



From the theoretical point of view the possibility of for- 
mation of local moments by increasing disorder has been 
suggested since the early developments of the theory of the 
combined effects of disorder and interaction [Q. In 3d the 
renormalization group (RG) flows to a Fermi liquid if dis- 
order is sufficiently weak. However, by suitably tuning the 
couplings of the RG one finds that indeed the combined 
effect of disorder and interaction leads to a magnetic in- 
stability signalling the formation of magnetized regions 0. 
An alternative point of view associates the local moments 
in 3d Si:P to the presence of rare fluctuations [|To|] . The 
relation and interplay between the two proposed mecha- 
nisms of local moment formation is still an open issue. 

Within the scaling theory the magnetic instability oc- 
curs, in 2d, for arbitrarily weak disorder and one argues 
that local moments form on a finite length scale at which 
the magnetic susceptibility diverges and the RG stops 0. 
2d systems then, apart from being presently of great in- 
terest, are good candidates for testing the RG approach 
to the interaction and disorder effects in the thermody- 
namics. 

In this paper we revisit the thermodynamic behavior of 
a two-dimensional disordered interacting electron system. 
In particular we show that by taking into account the cor- 
rection to the thermodynamic potential due to the com- 
bined effect of disorder and interaction in the particle-hole 
spin channel (triplet channel), there arises a non mono- 
tonic behavior of the specific heat both in temperature and 
magnetic field, strongly resembling the Schottky anomaly 
of free local moments. We interpret this as a dynamical 
precursor effect of formation of local moments. The above 
non monotonic behavior may well be difficult to observe 
by a direct measurement of the specific heat of the 2d sys- 
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tern. We then propose to perform magnetoconductance 
measurements in a time varying mag netic field and look 
for hysteresis (more details below) Q . We shall suggest to 
carry out this experiment well inside the metallic phase, 
for which we shall give predictions for the main energy 
scales. 

The starting point is the expression of the thermo- 
dynamic potential F in the presence of a magnetic field. 
We consider the corrections coming from the singlet and 
triplet particle-hole channels in the weakly localized regime, 
i.e. in lowest order in the dimensionless resistance. We ne- 
glect corrections from the particle-particle (Cooper) chan- 
nel and orbital effects The magnetic field is coupled 
via the Zeeman splitting of the spin states. These cor- 
rections were calculated in Ref.|lj| for weak interactions 
and extended to strong scattering amplitudes in Ref. jw| . 
Here, for completeness, we further extend these calcula- 
tions to include the energy renormalization Z which plays 
the role of m* jm in the context of the Fermi liquid the- 
ory of disordered systems Q , Q . The singlet and triplet 
particle-hole contribution to F is SF = J2j m F M wnere 
J and M are the total and the z component of the spin of 
the particle-hole pair (z being the direction of the applied 
magnetic field). J takes the values J = 0, 1 for the singlet 
and triplet contributions, respectively. 

In 2d the quantities F^j are 

r" 1 

F ° + F 1 = gN QP (-l + 72 ) Jdu (b (~) + ~) ^log (wt) 
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Here g = \/(Ait 2 NqD) — e 2 / '(wh)Rn is the dimensionless 
resistance in two dimensions, D is the diffusion coefficient, 
Nqp = ZN is the quasiparticle density of states per spin 
(N being the bare density of states), r is the elastic scat- 
tering time and b(u) — [exp (uj) — is the Bose func- 
tion. Eqs.(|-|) are valid in the case of long range Coulomb 
forces with 72 being the interaction coupling constant in 
the triplet particle- hole channel Finally, Q s — gLfJ-BH 
and Q s = (1 + 72)^ are the bare and interaction dressed 
Zeeman spin splitting frequencies, which enter the diffu- 
sive particle-hole propagators 

One can rewrite SF by separating the temperature and 
the magnetic field dependent parts as follows 



SF = F (T) + Fi(fl.) + gN QP T 2 f 2 



where 



F (T) = -gN QP (1 - 3 72 ) T [ — log (tT) + 



(3) 



(4) 



f 1 (n s ) = ± 9 n qp12 (i + l2 )n 2 io g n S T, 



(5) 



- 2 72 a (6) 



with a = Jp 00 dyyb(y) logy ~ —0.24 and 



f(x)= \ dyb(y) [(y - x) log \y - x\ + (y + x) log \y + x\] 
Jo 

(7) 

At small x, f 2 (x) k, -(1/2) 72 (1 + 7 2 )x 2 logs. 

At zero magnetic field, Eq.(Q) gives the leading loga- 
rithmic correction to the specific heat cy = —Td 2 F/dT 2 , 
which then is logarithmically enhanced at low tempera- 
tures for 72 > 1/3. This correction signals the break- 
down of perturbation theory for the quasiparticle den- 
sity of states and is related to the scaling equation for 
the energy renormalization Z || , |Q . The analysis of the 
RG equations in two dimensions shows that upon scal- 
ing both Z and 72 grow so rapidly that the renormal- 
ization procedure must stop at a finite length scale, L c = 
I exp[c/<7o(72,o + l)], where they both diverge. I is the mean 
free path, c a number of order one, and go, 72j o are the bare 
values of the running couplings g and 72 . Because Z and 72 
are related to the specific heat and to the spin susceptibil- 
ity via the relations cy/cy — Z and x/xo = ^(1+72), this 
divergence signals a magnetic instability and the above 
length scale L c has been interpreted as the typical size 
over which local moments form in the disordered system 
[|. The very existence of local moments is expected to 
have specific consequences in the behavior of thermody- 
namic quantities, in particular in the presence of a mag- 
netic field. In the case of free local moments, the specific 
heat shows the so-called Schottky anomaly, which man- 
ifests as a non monotonic behavior both in temperature 
and in magnetic field. On the experimental side, the be- 
havior in the specific heat and magnetic susceptibility has 
provided some support to the idea of formation of local 
moments in Si:P materials^. 

In the theory of disordered systems, a finite magnetic 
field enters as a mass term in the diffusive particle-hole 
triplet propagators with M = ±1, thus effectively cutting 
off the logarithmic singularities. As a consequence, the 
weak localization correction to the specific heat, Scy, is 
a decreasing function of the magnetic field H, at least 
when H 3> T. However, explicit evaluation of the leading 
(H/T) 2 term at small H/T shows that 5c v increases with 
the field, resembling the Schottky anomaly for free spins. 
It is useful to define the specific heat relative to the case 
with zero magnetic field 



Ac v (T, f2.) = cy(T, Q s ) - c v (T, 0) = -gN QP T 
I Ui[^- )--2~f' 2 [^ ) - ( ) 2 /2'(^M )(*) 




where f 2 and its derivatives can be evaluated numerically 
from Eqs.(^|-0). The result for Acy normalized to cyg^ 2 
with cy — (2/3)it 2 Nq P T is shown in Fig.l. 
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Fig. 1. Magnetic field dependent specific heat, 
Ac v (T, n s )/[c v 372] vs fl s /T. The full, dashed and long- 
dashed lines represent 72 = 1, 5, 10. 
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Fig. 2. Temperature dependence of the specific heat in pres- 
ence of a magnetic field Acv(T, f2 a )/[cvgj2](T/f2 s ) vs T/fi a . 
The full, dashed, long-dashed lines represent 72 = 1,5, 10. 

The quantity Acy(T, Q s ) has a non monotonic behav- 
ior as function of the magnetic field, resembling the Schot- 
tky anomaly of free local moments. The values f2 s ,mnx and 
fl s ,o at which Acv is maximum and zero, respectively, de- 
pend on 72 and move to smaller values upon increasing 72. 
A non monotonic behavior is also present in the tempera- 
ture dependence. This is shown in Fig. 2 where we report 
(Ac v (T, Q s )/(c v gj 2 ))(T/f2 s ) versus T/Q s . 

Introduction of a magnetic field results in Acv > 
(Ac v < ) for T > i? S:0 ( T < i? S)0 ). This should 
be contrasted with the case of free local moments where 
Acy > always. A magnetic field dependent specific 
heat is a direct consequence of the magnetic field de- 
pendence of the entropy. For free local moments this de- 
pendence leads to a non monotonic (dS/dH)T decreas- 
ing linearly at small H . We find that the excess entropy 
AS = S{T 7 H) - S(T,Q) which follows from Eq.(|) has 
a field dependence which mimics that of local moments. 
In Fig.3 we report dS(T,H)/dH vs H/T normalized to 
CV972/T. 

We also report dS\ oc (T, H)/dH for local moments. We 
rescale the density of local moments to fix the slope at 
small H to the slope of dS(T,H)/dH for 72 = 3. Like 
in the free local moment system, dS(T, H)/dH shows a 
minimum, whose location J7 Sjm i n depends on the strength 
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Fig. 3. Magnetic entropy (dS(T, n s )/dn a )T/ [cygn] vs O a /T. 
The full, dashed and long-dashed lines are numerical results for 
72 = 1,3, 10. The dashed-dotted line is the magnetic entropy 
for free local moments, where the density of local moments has 
been rescaled in order to fix the slope at small magnetic field. 

of the interaction. Notice that for 72 ~ 3 one obtains 
fts,min/T « 1.5, close to the value for free local moments. 

If local moments are thermally well coupled to the 
conduction electrons and weakly coupled to a heat bath, 
~(dS / 8H)t will control the rate of adiabatic "heating" 
and "cooling" of the conduction electrons in a time vary- 
ing magnetic field. The temperature hysteresis can then 
be made visible, e.g., by measuring the magnetoresistance 
p(H(t),T(t)). The change of temperature of the electrons 
is T-T = -(l/k)T{dS/dH) T (dH/dt)-(c v (T, H)/k)(dT/dt), 
where T is the temperature of the heat bath and k is 
a parameter describing the thermal coupling of the elec- 
trons with the bath. By using Ac v (T, H) and AS(T, H), 
one may estimate the order of magnitude of the hysteretic 
effect. Considering parameters which are in the domain 
of validity of the theory and correspond to the metallic 
phase, we take g = 0.05, 72 = 3, 10 and Z ~ 1. With a 
sweep rate of O.lTesla/minute and a thermal couling of the 
order k/cy ~ l/100sec we obtain a temperature variation 
of few (5-15%) percent indicating that the effect we are 
describing should be visible under the above parameter 
conditions. 

Concerning the size of the effects discussed in this pa- 
per and the possibility to observe them in a broader pa- 
rameter range one should keep in mind the following: 

i) Equations (0) and (||) have been derived by pertur- 
bation theory to lowest order in g (metallic limit). Indeed 
we can extend the region of validity of these expressions 
(and of Acy and AS) by assuming that g, 72 and Z are the 
running couplings of the RG analysis. In this case these 
couplings become temperature and field dependent and 
can assume large values, in particular Z> 1 even though 
572 < 1, which is the "optimistic" limit of confidence of 
the RG analysis 0,0] ■ This means that the predicted size 
of the hysteretic effect can be larger than the above con- 
servative estimates. 

ii) For f2 s < T the running couplings should renor- 
malize according to the RG equations derived at H = 
and are therefore only functions of the temperature. This 
implies that a measurement of Acy and of AS at small H 
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is a measurement of 72 (T) and Z(T). In particular we get 
Acy = -AS = {l/2)gN n TZ l2 {l + l2 )(n s /T) 2 at leading 
order in (H/T) 2 . It would be useful to measure these quan- 
tities in the metallic side of 2d systems showing a MIT to 
assess the validity of the interacting scaling theory of the 
2d metallic behavior B. 

iii) The occurrence of a local moment instability at 
finite length makes relevant the sample dishomogeneity. 
Small local variations of disorder could result into regions 
with local moments coexisting with "more" metallic re- 
gions. Here, with regions of local moments we mean re- 
gions where the interplay of disorder and interaction in 
the spin channel has reached the strong coupling limit 
gi2 > l 7 Z ^> l,x/Xo ^ 1- On a physical ground, one 
indeed expects that the weak coupling RG description of 
local moments in the metallic phase will evolve at low tem- 
perature into a more localized picture similar to that of 
the insulating phase |19[| , possibly in terms of an effective 
Kondo Hamiltonian pOt ^ However a comprehensive theory 
of this crossover from weak to strong coupling is still lack- 
ing and is a main open problem. A systematic analysis of 
the validity and of the breakdown of the existing scaling 
theory of disorder and interaction can shed light on this 
issue. 

We like finally to comment on 3d systems. The 2d 
results in Eqs.(Q) and (||) are easily extended to d > 
2. Acy and AS show non monotonic behavior as func- 
tion of T and fi s analogous to that in 2d. The leading 
(H/T) 2 contribution to Ac v and AS now reads Acy = 
0.23gZ^ 2 (l + l2 )^ 2 (^T+^-l)N T^/T/{K^(as/T) 2 . 
and AS = —2Acy. By assuming that this expression for 
Acy also holds approaching the strong coupling regime, 
where both x/xo and cy/c v diverge, one obtains {Acy / 
{Q s /T) 2 ) ~ (x/xo) 2 (c Q y/cy)^ 2 T 3 / 2 . Also this result can 
be tested experimentally to fix the limits of validity of the 
weak coupling description of local moments in 3d systems. 

C.C. thanks M.Fabrizio for various helpful discussions. 
This work was supported by the TMR program of the 
European Union (P.S.). 
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